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Antiferromagnetic (AFM) materials are widely used in spintronic devices as passive elements 
(for stabilization of ferromangetic layers) and as active elements (for information coding). In both 
cases switching between the different AFM states depends in a great extent from the environmental 
noise. In the present paper we derive the stochastic Langevin equations for an AFM vector and 
corresponding Fokker-Planck equation for distribution function in the phase space of generalised 
coordinate and momentum. Thermal noise is modeled by a random delta-correlated magnetic field 
that interacts with the dynamic magnetisation of AFM particle. We analyse in details a particular 
case of the collinear compensated AFM in the presence of spin-polarised current. The energy dis- 
tribution function for normal modes in the vicinity of two equilibrium states (static and stationary) 
in sub- and super-critical regimes is found. It is shown that the noise-induced dynamics of AFM 
vector has peculiarities compared to that of magnetisation vector in ferromagnets. 
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I. INTRODUCTION 

Magnetic nanoparticles are the main constituents of the nowadays devices for information technology. While the 
deterministic dynamics of magnetisation vectors is used for information coding, the noise- induced stochastic behaviour 
facilitates the switching processes and thus is used to increase the speed of information processing (see, e.g.^'"). Noise 
measurement is a powerful and informative tool for study of the spintronic effects in different systems '" '^. The details of 
the stochastic processes are also important for the development of high-quality spin-torque oscillators ' and micropower 
generators ' . 

Theoretical approach to the description of thermal noise in small ferromagnetic (FM) particles was developed in the 
seminal papers of W.F. Brown ' ' where the thermal bath was modeled by fluctuating magnetic field and corresponding 
Langevin equations were obtained as generalisation of the dynamic Landau-Lifshitz-Gilbert equations. This approach 
was then extended to the systems in the vicinity of the Curie point , space- inhomogeneous magnetic vortices^ ^, 
systems with the coloured noise^''^ ' and FMs in the presence of spin-polarised current'"''' '. 

On the other hand, stochastic behaviour of antiferromagnetic (AFM) nanoparticles which are also widely used in 
electronic and spintronic devices is studied to a much lesser extent. The main reason is in the seeming magnetic 
neutrality of AFMs which manifests itself in the vanishingly small or zero net magnetisation, quadratic (in contrast 
to linear for FMs) dependence of the internal energy vs external magnetic field etc. Thus, up to now the problem 
of magnetic relaxation and thermal noise in AFM particles was in fact reduced to the description of the effective 
FM particle with small but nonzero magnetisation which inevitably appears due to imperfections, strong external 
magnetic fields, surface effects etc"'"' '. However, in many cases a peculiar feature of AFM, namely, the presence 
of strong exchange coupling between the differently aligned (mainly in opposite directions) magnetic moments, gives 
rise to new dynamic effects that could not be reduced to the motion of the above mentioned "stray" magnetisation. 
For example, in contrast to magnetisation of FM particle, an AFM (Neel) vector (defined as the difference between 
sublattice magnetisations) can be set into motion not only by external magnetic field H but also by its time derivative 
H (or curled electric field) . Typical frequencies of AFM oscillations are fall into terahertz range (compared to 1-10 
GHz for FM) due to effects of the exchange enhancement. AFM systems, like FMs, are sensitive to spin-torques 
transferred from the spin-polarised current'^'"- ' but the current-induced AFM dynamics differs significantly from the 
current-induced dynamics of FMs- '. Thus, with the account of perspective to use AFM nanoparticles as alternative 
to FM active elements of spin- valves (see, e.g. - '), theoretical study of thermal noise in these systems is of importance 
and of great interest. 

In the present paper we generalise the dynamic equations for AFM nanoparticle to the stochastic Langevin equations 
that describe "Brownian motion" of AFM vector in the presence of thermal noise. Like in the W. F. Brown approach 
we model the noise as a fluctuating delta-correlated magnetic field which interacts with the dynamic (induced by the 
motion of AFM vector-'"'- ' ) magnetisation of AFM. Using the standard Langevin dynamics technique for multiplicative 
noise-^ we derive the Fokker-Planck equations for the distribution function in the phase space of AFM vector and 
corresponding generalised momentum and discuss the peculiarities of AFM system compared to FM one. As a 
representative example we analyse the stochastic behaviour of AFM particle in the presence of spin-polarised current 
and find the energy distribution function in subcritical and supercritical regimes. For the sake of simplicity we restrict 
ourselves with the case of a collinear compensated AFM with two oppositely directed magnetic sublattices. However, 
the developed approach allows generalisation for the multisublatteral AFM, weak ferromagnets etc. 

II. LANGEVIN EQUATIONS FOR ANTIFERROMAGNETIC PARTICLE 

In what follows we consider the fine (nanosized) magnetic particles whose state can be described with only few 
macroscopic vectors M^: magnetisation vector (j = 1) in the case of FM nanoparticle and vectors of two equivalent 
sublattice magnetizations {j = 1, 2) in the case of a collinear AFM. 

Let us first discuss the stochastic Landau-Lifshits-Gilbert equation for the FMs which describes the dynamics of 
magnetisation M subjected to a spin transfer torque Tstt at a finite temperature T^'^: 

M = 7M X (Heff + h) + X [M X (H^ff -I- h)] + Tstt, (1) 

where 7 is the gyromagnetic ratio, Hofr = — 9i(Jan(M)/(3M is the effective (combination of internal and external) 
magnetic field, Wan(]V[) is the magnetic anisotropy energy with Zeeman contribution, ac is the dimensionless damping 
(Gilbert) constant, the sign x means cross-product. The fluctuating magnetic field h(i) with a Gaussian stochastic 
process has the following standard space-time statistical properties: 



(h(i)) = 0, {h,{h)hk{t2)} = 2DSjkS{h - t2), 



(2) 
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where D represents the strength of thermal fluctuations whose value is defined from the fluctuation-dissipation theo- 
rem. Symbol (. . .) denotes an average taken over all realizations of fluctuating field. 

The first r.h.s. term in Eq. (1) describes nondissipative dynamics of magnetisation which is analogue of the precession 
of gyroscope with the "angular momentum" M. The second r.h.s. term describes dissipation processes that have the 
same origin as fluctuating field h.(t). Both nondissipative and dissipative terms are subjected to noise, however, the 
nondissipative noise is additive while the dissipative one is multiplicative . Delta-correlation of fluctuating field (2) 
means that the noise correlation time is much lesser than the characteristic time of magnetisation response. 

The dynamics and kinetics of FM magnetisation with account of normalisation condition |M| = Mq (imposed far 
below the Curie point) is described with two independent variables that define space orientation of vector M, in other 
words, by the variables of configuration space. Full "mechanical" energy of FM particle, EpM = u;an(M), also depends 
upon orientation of M and thus can be treated as consisting of the potential energy only, in contrast to the energy of 
AFM particle (see below). 

Using the analogy with FMs, one can derive the stochastic equation for AFMs from the corresponding dynamic 
equations assuming that the thermal noise also has the magnetic nature and can be modeled with the same random 
field h.{t) (2). This field may originate from fluctuations of i) the surface noncompensated magnetisation for the 
small particles; ii) magnetisation of the nearest FM layer in spin-valves and multilayers; in) current that produces 
additional magnetic field. 

However, the deterministic dynamics in AFMs substantially differs from that of FM magnetisation and looks like 
an inertial motion of a point mass in a potential well. Formally this effect was demonstrated for AFMs with strong 
exchange coupling between the magnetic sublattices"'*'"' . In this case the dynamics of a collinear AFM is described 
by a single AFM (so called Neel) vector L = Mi - M2 of a fixed length (|L| = 2Mq). 

Deterministic equations of motion for L could be obtained either from the set of Landau-Lifshits-Gilbert equations 
(I) for each of magnetic sublattices or, equivalently, in the framework of Lagrange formalism (see ' for details). The 
last approach is more convenient for general analysis, so, we start from the Lagrange function for AFM particle in 
the following form 



jOafm = — L + -yniL 



L.(LxH) -i«,„(L) + l^(LxH)2, (3) 



where H is an external magnetic field and Wan(L) is the energy of magnetic anisotropy that forms a potential well for 
AFM vector, 7, as above, is the gyromagnetic ratio. The value = 1/{2j'^MqHe) plays a role of the "inertia mass"; 
it depends upon the spin-flip field of the exchange nature, 2He, that characterises the intersublattice coupling. 

In contrast to FMs, the nondissipative dynamics of AFMs can be also described within the Hamiltonian formalism 
with the following generalised energy (Hamilton function ) obtained from (3): 

Hafm = TT^Pi - 7 [Pi ■ (L X H)] + Wa„(L). (4) 

Here the generalised momentum Pl = 9£afm/9L is canonically conjugated to the generalised coordinate L. 

As can be directly seen from (4), the generalised energy of AFM particle, Hafm, includes both kinetic (first term) 
and potential (last term) contributions. It means that the dynamics and kinetics of this system is described within 
the phase (vs configuration for FMs) space. 

Dissipation is modeled with the Raileigh function which in the presence of spin-polarised current J takes a form-"': 



-T, f 2 

K-AFM — 7AFMnT-L-L — „ , . 

27M0 



(LxL) . (5) 



Here the first term models the internal damping, damping coefficient 27afm is the AFMR linewidth, the constant 
a — h^e / {2eMQVAFM) is proportional to the efficiency e of the spin transfer processes, uafm is the volume of AFM 
nanoparticle, h is the Plank constant, e is the electron charge. Unit vector Pcuir is parallel to the direction of the 
current spin polarisation. 

Thus, the stochastic equations for AFM in the phase space {L,Pl} obtained from (3) and (5) with substitution 
H — > h{t) acquire the form: 

L = Pl/tol - 7L x h (6) 

Pl = Fl + Fdiss - 7 (Pi - 27afm'71lL) x h, 

where = — (9i(;an(L)/9L is the potential (gradient) force, and the dissipative force Fdiss is given by the following 
expression 



^ _ 97?.AFM 
diss — 



dL 



= -27afmPl - Pcurr X L. (7) 
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Equations (6) describe the evolution of AFM vector in the presence of thermal noise and in this sense are analogous 
to the stochastic Eqs. (1) for magnetisation M of FMs. Both sets of equations (for both FMs and AFMs) are linear 
in the random magnetic field h. For an AFM system this fact is nonobvious and can be explained by the presence 
of small but nonzero macroscopic magnetisation Mafm = Mi + M2 that in the compensated AFM has a dynamic 
origin and can be expressed in terms of the Neel vector: Mafm ocLxLocPi xL. On the other hand, the noise 
terms in both (FMs and AFMs) equations are multiplicative and this can, in principle, result in a possible stochastic 
resonance. 

It should be stressed that though the AFM dynamics is similar to the dynamics of point mass, Eqs. (6) have one 
peculiarity compared with the standard Langevin equations for a Brownian particle in a potential well. Namely, the 
first of Eq. (6) includes the noise and does not include any dissipation term. This means that within the accepted 
model of dissipation (and noise) there is no time-scale separation between relaxation of generalised coordinate and 
generalised momentum. This fact is a direct consequence of limitations on |L| imposed by assumption of strong 
exchange coupling and absence of the exchange relaxation. However, the characteristic energy of exchange coupling 
is of the order of the Neel temperature. So, for low (compared with the Neel) temperatures and relatively small 
(compared with He) external fields Eqs. (6) give an adequate description of AFM vector behaviour. 

The Langevin Eqs. (6) generate the Fokker-Planck equation for AFM probability distribution function /(L,Pi;t) 
in the phase space: 



+ Vp, • 



2MoHePl + Dk^ • Vl - Pl ® L ■ Vp^ j / 
Fl + Fdiss + D-i'^k^^ • Vp^ - L ® Pl • Vl) / , 



(8) 



where we introduced the symbol A"* = la? — a ® a with a = L or P^, and omitted small noise terms with 7afm for 
the sake of clarity. 

Fokker-Planck Eq. (8) for AFM nanoparticles is, in fact, the main result of this paper. It is much more complicated 
than the analogous equations for FMs and (in contrast to FMs) could not be solved in general case even for stationary 
conditions. In the next section we consider some limiting cases that allow to find approximate stationary solutions, 
/(L, Pi), and evaluate D from fluctuation-dissipation theorem. 



III. ANTIFERROMAGNET PROBABILITY DISTRIBUTION IN THE PRESENCE OF 

SPIN-POLARISED CURRENT 



The phase space of AFM particle is four-dimensional and this substantially complicates analysis of Eq. (8) in 
general case. However, in some cases the effective dimensionality can be reduced. The simplest case concerns the 
system in the vicinity of equilibrium where all the possible motions of AFM vector could be represented in terms of two 
noninteracting normal modes with the amplitudes c± and angular phases Lp± . If, in addition, we neglect inhomogenuity 
in the phase Lp± distribution, then, distribution function can be factorized as: /(L, P^,; t) = /_|_(c_|_; t)f-{c-] t). In what 
follows we consider the case of AFM with the degenerate excitation spectra for which two normal modes correspond 
to clockwise/counter-clockwise rotations of AFM vector around z axis with the frequency r^AFMR (that is close to 
AFMR frequency). 

Spin current polarised along z axis (pcurrll-z) interacts with both modes thus enhancing the effective damping of 
one (say, "-I-") and diminishing the effective damping of the other (say, "-")- '. In the subcritical regime (|J| < 
Jcrit = '^Iafm^afmrI (n<^ He)-, positive damping) the static equilibrium state is stable and normal modes are still 
well separated. In the supercritical regime, \J\ > Jcrit, an amplitude of one of the mode growth to saturation 
value and the stable state corresponds to rotation of AFM vector in xy plane with the current-dependent frequency 
Lu = Jf^AFMR/'^rit'"'- Another normal mode corresponds to small oscillations of AFM vector in z direction, so, again, 
both modes are well separated. Thus, the behaviour of AFM vector in the subcritical and supercritical regions can 
be really described in approximation of two independent normal modes. 

To obtain the Fokker-Planck equations for f{c±) we use the approach of energy representation for nonequilibrium 
Brownian- like systems developed in". To this end let us start from the Langevin equation for the energy Eafm = 
P|/(2mL) + Wan(L) (compare with (4)): 



dt 



Pr -Fh 



+ 7 



27. 



di 



afmJ 



L X h 



(9) 



where the summands in the r.h.s. of Eq. (9) should be expressed in terms of Eafm- 
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In approximation of noninteractive normal modes Eq. (9) is applicable to the energy E± = 4Mq fi^pj^j^miCj. of 
each mode. Moreover, within the accepted approximation (fixed oscillation frequency riAFMn) E± could be considered 
as the dynamic variables (that are proportional to the "true" canonical variables, actions). 

In the subcritical region, \J\ < Jcrit, Eq- (9) can be rewritten as follows: 

—J- = -7AFM 1 ± -J C± - 7- Ci,2/lz (10) 

V Jcrit/ "AFMR 

+ 27- {ha; cos r^AFMRi + hy sin r^AFMRi) 

"AFMR, 

+ 7 {hoc sinf^AFMRi - hy cosOafmrO . 

It should be stressed that the same equation could be obtained directly from (6) after transition to amplitude-phase 
representation. 

As it is seen from Eq. (10), the sign ± corresponds to different modes which interact with the current in different 
ways. If J > 0, the effective damping of the first mode (with the amplitude c+) increases and that of the second (with 
the amplitude c_) decreases, due to the action of spin-polarised current. 

Analysis of Eq. (10) shows that one component of the random magnetic field h, namely, that, which is perpendicular 
to the plane of L rotation {hz in our notations), is a source of multiplicative noise. However, if the damping is rather 
small, 7AFM ^^AFMR, the term with multiplicative noise can be omitted. To this end Eq. (10) generates the 
following Fokker-Planck equations: 



df{c±) ^ d 
dt dc± 



7AFM 1 ± -J — c±/(c±) + Dj 



Jrrit, / dC4 



(11) 



From the stationary solution of (II) one gets the AFM probability distribution function f{E+,E-): 

(12) 



f{E+,E_) = /oCxp<^ -- 



^^AFMR^-^O 



i + ^]e^+(i-J-\e_ 



where /o is a normalization constant. 

In the absence of current the distribution (12) should coincide with the Boltzmann distribution function. From the 
fluctuation-dissipation theorem we get the diffusion coefficient for AFM particle 

„ IAFmHe rr 1 7AFM I^E rj. (^n\ 

^ ^AFMR-'^'-'O 7^"0 "AFMR V J^a 

where Ha = ^^afmr/(7^-^s) is the field of magnetic anisotropy which in the typical AFMs is small compared with 
strong exchange field: Ha ^ He- 

Remind, that the analogous coefficient for FM particle has a form'^"': 

Dfm = ^^T, (14) 

where we used an explicit expression for the Gilbert damping parameter through the frequency and half-width of 
FMR, aa = 7FM/i^FMR- Comparing (13) and (14) one can easily see that the diffusion coefficient in AFMs is greater 
that that for FMs due to the large factor HE/Ha ^ 1 , other things being equal. This is one more manifestation of 
the above mentioned exchange enhancement peculiar to AFM materials. 

In the presence of spin-polarised current the distributions (12) are still Boltzmann-like with two (instead of one for 
FMs) different effective temperatures for each mode: 

1 ± J/ Jcrit 

Expression (1-5) shows that the temperature of the "soft" mode (that one which becomes unstable at J ^ Jcrit) 
crucially growth, while the temperature of the other mode diminishes. This fact illustrates the current-induced 
energy swap between two modes. Seeming singularity at J — Jcrit is an artifact of approximation which presupposes 
existence of high energy barrier between the different stable states. 

In the supercritical region one can get the distribution function in a similar way. Neglecting, whenever it is possible, 
the small value 7AFM/f^AFMR ^ 1, we arrive at the following expression: 



/(S+,i;_) -/oexp 



AE^ 



E^ - 



[3 + 4( J/ Jerit)2]T 2rii;i") 



(16) 
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where is related with oscillations of AFM vector in z direction and the average energy of the second mode (related 
with the rotation of AFM vector in xy plane), = M^HaiJ / Jcntf is proportional to the current value. 

Like in FMs'', the distribution (16) is Gaussian-like with respect to the energy of the second mode. However, in 
contrast to FM, the half-width of corresponding distribution is proportional to Ai?_ oc J, so the "quality factor" 

e'^'/AE- oc J growth with the current value . Another peculiarity of AFM system compared with FM is the 
presence of the additional energy fluctuations related with the first mode. 



IV. CONCLUSIONS 

In the present paper we have derived the Langevin and the Fokker-Planck equations that take into account the 
peculiarities of the dynamics of AFM (in contrast to magnetisation) vector. These equations could be used for 
calculations of the dwell times between the different states of AFM particle and the lineswidth of resonances induced 
by external fields (including spin-polarised current). It is shown that the thermal noise generated by fluctuating 
magnetic field is multiplicative. As a result, corresponding Fokker-Planck equation is nonlinear and this opens a 
possibility for noise-induced transitions and stochastic resonances in the system. 

In the framework of the proposed approach we have calculated the AFM energy distribution function for the 
particular case of the collinear two sublattice AFM in the presence of spin-polarised current. We found that at a 
given temperature and quality factor of the magnetic resonance the diffusion coefficient of AFM shows an exchange 
enhancement compared to that of FM nanoparticle. It is also shown that spin-polarised current affects the effective 
temperatures of the normal oscillation modes in different ways: in the subcritical region the temperature of the soft 
mode increases and the temperature of the other ("hard") mode decreases. In the supercritical region the energy 
fluctuations of the soft mode grow with respect to the current value slower than the average energy of the mode. This 
opens a way to control the efficiency of energy transfer from the current to AFM oscillator. 

In our modeling we considered only the magnetic sources of noise. However, in the presence of spin-polarised current 
the fluctuations of the current value could be a source of multiplicative noise, as seen e.g. from Eq. (10). The problem 
of current-induced noise needs a special treatment that accounts for the relations between the magnetic state of AFM 
layer and resistivity. Joule losses etc. 

Another important extension of the problems considered in the present paper is seen in analysis of the possible 
current-induced nonequilibrium states and their thermodynamics and information characteristics in the spirit of recent 
general approaches"^'" ' ' . 
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